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Abstract: Forest recovery may be influenced by several factors, of which fire is the most critical. However, 
moderate- and long-term effects of fire on forest recovery are less researched in Northwest China. Thus, 
the effects of different forest recovery time after fire (1917 (served as the control), 1974, 1983 and 1995) 
and fire severities (low, moderate and high) on larch (Larix sibirica Ledeb.) forest were investigated in the 
Kanas National Nature Reserve (KNNR), Northwest China in 2017. This paper analyzed post-fire 
changes in stand density, total basal area (TBA), litter mass, soil organic carbon (SOC) and soil nutrients 
(total nitrogen, total phosphorus and total potassium) with one-way analyses of variance. Results indicate 
that litter mass, TBA, SOC and soil nutrients increased with increasing recovery time after fire and 
decreasing fire severity, while the stand density showed an opposite response. The effects of fire 
disturbance on SOC and soil nutrients decreased with increasing soil depth. Moreover, we found that the 
time of more than 43 a is needed to recover the litter mass, TBA, SOC and soil nutrients to the pre-fire 
level. In conclusion, high-severity fire caused the greatest variations in stand structure and soil of larch 
forest, and low-severity fire was more advantageous for post-fire forest stand structure and soil recovery 
in the KNNR. Therefore, low-severity fire can be an efficient management mean through reducing the 
accumulation of forest floor fuel of post-fire forests in the KNNR, Northwest China. 
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1 Introduction 


Wildfire is an important disturbance agent in forest ecosystems (Liu et al., 2013; Alcaniz et al., 
2018), especially in boreal forests (Simard et al., 2001). Various studies have focused on the fire 
regime, an important component of which is fire severity (Thode et al., 2011). Fire severity can be 
defined as the total effect of a fire on an ecosystem (White and Pickett, 1985; Savadogo et al., 
2007), including tree damage and mortality as well as the quantity of burnt organic matter (Wright 
and Clarke, 2007; Jayen et al., 2016). Fire severity affects soil properties (Certini, 2005), post-fire 
vegetation community structure and dynamics (Certini, 2005; Wang and Kemball, 2005; Cai et al., 
2013; Jayen et al., 2016), competitive interactions among plants in post-fire environments (Ducey 
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et al., 1996), post-fire recovery levels in numerous forest systems (Perevoznikova et al., 2007; 
Vivian et al., 2008; Knelman et al., 2015), and it can also change plant species diversity and 
richness through changes in site characteristics (Aref et al., 2011; Pourreza and Hosseini, 2014). 
Low- and moderate-severity fires increase the diversity and richness of plant species (Pourreza 
and Hosseini, 2014); however, high-severity fire has been shown to exert an opposite effect 
(Pywell, 2002; Garcfa-Dominguez and Fernandez-Palacios, 2009). For example, the number of 
indicator species is known to decrease with increasing fire severity (Johnson, 1992; Heydari et al., 
2016). 

Fire led to large changes in stand structure (Peterson and Reich, 2001), and physical, chemical 
and biological properties of soils (Gonzdlez-Pérez et al., 2004; Certini, 2005). Burned sites 
generally had lower tree density than unburned sites (Signell et al., 2005). Fire causes increases in 
soil pH, bulk density and total phosphorus (TP) (Fernandez et al., 1997; Martin et al., 2012) and 
decreases in soil organic matter and carbon (C) (Neary et al., 1999; Johnson and Curtis, 2001; 
Certini, 2005), thus altering the soil physical structure (Arocena and Opio, 2003). 

The post-fire recovery of plant species diversity needs a long-term process (about 90-100 a) 
(Zyryanova et al., 2007). Low-severity fire is advantageous to forest development (Bergeron, 2010). 
Post-fire recovery of soil C, N (nitrogen) and P (phosphorus) relies on organic matter 
accumulation that linked with forest productivity (Quintero-Gradilla et al., 2015). Furthermore, 
forest fire severity impacts recovery rate. The recovery rate of forest following a low-severity fire 
is higher than those of following moderate- and high-severity fires. Concentrations of C and N in 
Nothofagus pumilio Krasser after 4 a post-fire were 52% lower than those of without fire (Alauzis 
et al., 2004). It may take 80 a for forest litter to recover to the pre-fire level (Lorenza et al., 2000). 
Understanding the post-fire responses of different soil layers is extremely essential for forest fire 
management. However, most studies have focused on the short- or moderate-term effects of fire 
severity (Gould et al., 2002; Lezberg et al., 2008; Kishchuk et al., 2015; Francos et al., 2016). A 
long-term perspective of post-fire responses should enable improved understanding of the 
resilience and capacity of a forest to respond to fire disturbance (Francos et al., 2018). 

The immediate effects of fire on forest have been widely documented (Quintero-Gradilla et al., 
2015), thus, the moderate- and long-term effects of different fire severities should be further 
investigated. Little is known about stand structure and soil property responses to fire severity in 
the Kanas National Nature Reserve (KNNR), Northwest China. Thus, the objective of this study 
was to analyze the changes of stand structure, litter mass, SOC, total N (TN), total P (TP) and 
total potassium (TK) in different soil layers after three times of fire occurrence following low-, 
moderate-, and high-severity in the KNNR, Northwest China. It was hypothesized that (1) the 
decreases of total basal area, litter mass, SOC and soil nutrients (TN, TP and TK) in high-severity 
fire stands were more significant than those of moderate- and low-severity fire stands; and (2) the 
recovery rates in low-severity fire stands would be the fastest among the three times of fire. 


2 Materials and methods 


2.1 Study area 


The study was carried out in the KNNR (48°36'18"—48°38'56"N, 87°01'45"—87°33'50"E; a total 
area of 875 km), located in Buerjin County, Xinjiang, Northwest China. The KNNR has a typical 
temperate climate. The annual average temperature is 0.2°C, the average annual precipitation is 
1065 mm, the average annual evaporation is 1097 mm, and the frost-free period is approximately 
80 to 108 d (Liu et al., 2009). 

There are 83 families, 298 genera and 798 species of vascular plants in the KNNR; among 
these, the dominant tree species are Larix sibirica Ledeb., Picea obovata Ledeb., Pinus sibirica 
(Loud.) Mayr., Abies sibirica Ledeb., Betula pendula Roth. and Populus tremula Linn. (Liu et al., 
2009). The dominant shrub species are Spiraea chamaedryfolia L., Rosa spinosissima L. and 
Lonicera cearulea L. The dominant herbaceous species are Carex atrata Linn., Galium boreale 
Linn., Aconitum leucostomum Worosch. and Elymus sibiricus L. (Fang et al., 2014). The study 
area has a long history of fire, most of which are originated from lightning. 
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2.2 Sampling design 


Field surveys and soil sampling were conducted from June to September 2017. The larch (L. 
sibirica) that is the dominant tree species in the KNNR, thrives in post-fire habitats 
(Gonzalez-Perez et al., 2004; Kuz'Mina, 2004), and their fire scars provide important data for 
studying fire regimes. Age analysis of forest fire scars was conducted to determine the years of 
fire disturbance (Farris et al., 2013) and tree rings were used to reconstruct fire history. Fire 
severity was determined using the height of tree scorch (Knapp and Keeley, 2006; Lentile et al., 
2006; Halofsky and Hibbs, 2009; Keeley, 2009), and the height of tree scorch refers to the 
distance from the ground to the height of tree blackening, measured with steel tape or a tree 
height measuring instrument LD6172. 

Larch stands that had been disturbed by fires in 1974, 1983 and 1995 were investigated. Stands 
that had only trace sign of a fire in 1917 served as the control (unburned stands). Study plots were 
stratified by fire severity (low-severity (LS), moderate-severity (MS) and high-severity (HS)) and 
by the year of the fire occurrence (1974, 1983 and 1995). To compare the effects of fire on stand 
structure and SOC, TN, TP and TK in three soil layers with recovery time after fire and fire 
severities, we used a space-for-time substitution design in this study that relied on the assumption 
that all variations among stands were due to differences in time since the fire disturbance (Yanai 
et al., 2003; Bergeron, 2010). 

By setting 100-metre-wide transects along altitude, fire-scarred trees were identified, and plots 
were established in larch stands with at least five fire-scarred trees at the same time. To minimize 
the influence of site conditions, we chose plots within the same altitudinal ranges from 1685 to 
1807 m in the middle of sunny slopes (Table 1). We selected 93 plots based on accessibility. The 
plots were 30 mx30 m in dimensions and the boundaries of the plots were at least 50 m from the 
forest edge. 


Table 1 Main characteristics of the studied plots 


ov p Year of fire Recovery No. of : 5 ; Canopy density 
Fire severity occurence time (a) plots Altitude (m) Slope (°) Aspect (%) 
1974 43 10 1696-1796 iG. Northwest, 60-70 
Northeast 
HS 1983 34 11 1685-1748 icio. Northwest; 38-50 
Northeast 
1995 22 7 1685-1759 1-10 Northwest 57-10 
1974 43 13 1694-1807 jig- Northwest, 56-70 
Northeast 
MS 1983 34 9 1693-1803 1-10 Northwest 38-54 
1995 22 10 1690-1741 1-10 Northwest 65-70 
1974 43 7 1680-1796 1-10 Northwest 60-65 
LS 1983 34 9 1698-1806 1-10 Northwest 53-70 
1995 22 8 1688-1807 1-10 Northwest 54-70 
Control 1917 100 9 1688-1777 1-10 Northwest 54-76 


Note: LS, low-severity; MS, moderate-severity; HS, high-severity. 


In each plot, we measured the diameter at breast height and canopy density. The altitude, 
gradient, aspect and slope of the plots and the species and number of all woody plants were 
recorded. Then, the TBA and stand densities of each plot were calculated. Five small quadrats (1 
mxl m) were set in the centre and four corners of each plot. Litter was collected from each 
quadrat, weighed and dried, and the litter mass was calculated. 


2.3 Soil sample collection and measurement 


To obtain soil samples, we investigated three soil profiles on the upper, middle and lower slopes 


in each plot. One-kilogram soil samples were collected from three layers (0—10, 10-25 and 25-70 
cm) in each plot and then mixed to obtain composited samples for each soil layer (Marañón et al., 
1999). Soil bulk density was sampled using 100-cm? circular knives in different layers, and 
samples were then air dried and sieved prior to analysis. Soil pH was determined in aqueous soil 
extract (1:2, soil:water) using a digital pH metre. Soil organic material was assessed using the 
potassium dichromate oxidation heating method (Walkley and Black, 1934). TN concentration 
was determined with the semimicro-Kjeldahl method (Bremner and Mulvaney, 1982), TP 
concentration was determined by the Bray extraction method (Bray and Kurtz, 1945) with a 
molybdenum extraction colorimeter, and TK concentration was determined by ammonium acetate 
extraction (Knudsen and Petersen, 1986) with the flame photometric method. The concentrations 
of SOC, TN, TP and TK multiplied by soil depth and bulk density of their areas of fire are equal 
to their contents. 


2.4 Data analysis 


One-way analyses of variance was used to test the effects of recovery time after fire and fire 
severities on each variable, and post hoc Tukey's tests (P<0.05) were performed for comparisons 
between means. These statistical analyses were carried out with SPSS version 24.0. Linear 
regression was used to assess the relationships between SOC and TBA, SOC and litter mass, and 
litter mass and TBA using OriginPro 2018. The structural equation model (SEM) was applied to 
analyse the effects of fire disturbance on SOC and soil nutrients using SPSS AMOS 23.0. 


3 Results 


3.1 Stand structure, soil property and litter mass 


There were significant differences in both stand density and TBA among the plots under different 
recovery time after fire and fire severities, suggesting that the coupling of forest fire severity and 
recovery time after fire strongly altered the stand structure of larch forests (Table 2). TBA 
increased with increasing recovery time after fire and decreasing fire severity, while stand density 
showed opposite trends and the stand density of the control was not smaller than those of the LS 
plots in 1974. 


Table 2 Stand density and total basal area among the plots under different recovery time after fire and fire 
severities 


LS MS HS 
1917 


Variable (Control) 1974 1983 1995 1974 1983 1995 1974 1983 1995 
(43 a) (34 a) (22 a) (43 a) (34 a) (22 a) (43 a) (34 a) (22 a) 


Bland density 1264 1175 1237? 1338 1336% 1385 1394% 1428? 1446 1571 


(trees/hm?) 
SD (trees/hm?) 54 84 16 70 63 93 74 64 90 78 
Total Dasal arsa 47.0° 44.2% 43.9% 41.8? 34.1% 351% 34.9% 34.9% 30.7 25.2 
(m?*/hm7) 
SD (m?/hm?) 0.5 5.4 11.8 9.6 2.7 3.0 5.9 3.4 1.8 1.3 


Note: SD, standard error. Different lowercase letters within a row indicate significant differences among different recovery time after 
fire and fire severities at P<0.05 level. LS, low-severity; MS, moderate-severity; HS, high-severity. 

Recovery time after fire and fire severity had significant effects on soil physical and chemical 
properties (Table 3). There were no significant differences in bulk density and pH among the plots 
under different recovery time after fire and fire severities. The soil pH values were moderately 
acidic (5.32-5.91). There were significant differences in SOC and TN concentrations between the 
control and the plots with different recovery time after fire and fire severities. The control plot 
had significantly higher SOC concentration than the HS plot in 1995 and had significantly higher 
TN concentration than the MS and HS plots in 1983 and 1995 (Table 3). The control plot had the 
highest SOC, TN, TP and TK values. There were no significant differences in the average 
concentrations of SOC, TN, TP and TK among the plots under different recovery time after fire 
and fire severities (Table 3). There were significant differences between the C:N and N:P ratios. 
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The C:N ratio of the LS plots in 1974 was significantly lower than those of the MS plots in 1974 
and 1983 and the HS plots in 1995, and the N:P ratio of the LS plot in 1974 was significantly 
lower than that of the control (Table 3). 


Table 3 Soil physical and chemical properties among the plots under different recovery time after fire and fire 


severities 


ine LS MS HS 
Parameter control) 1974 1983 1995 1974 1983 1995 1974 1983 1995 
(43 a) (34 a) (22 a) (43 a) (34 a) (22 a) (43 a) (34 a) (22 a) 
Bulk density 1991" 1.159" 1.062" 1.026° «1.033 1.093? 1.096° 1.006 1.091° 1.066 
(g/100 cm”) 
g e 0.466 0.064 0.068 0.093 0.062 0.029 0.069 0.987 0.506 0.055 
pH 5.290" 5.770"  5.510° 5.280 5.280 5.670 5.300° 5.700° 5.880" 5.200? 
SD 0.080 0.170 0170 0.120 0.160 0.180 0.100 0.150 0.250 0.130 
SOC (g/kg) 25.365" 19.340% 15.559% 14.126" 17.708 15.563% 14.126% 15.590% 15.138 12.147? 
SD (g/kg) 4.571 4.056 2.482 1.527 1.456 2.498 1.527 1.616 2.184 1.173 
TN (g/kg) 1.512" 1.244 0.934% 0.973% 1.011% 0.872? 0.819 0.911% 0.862? 0.668? 
SD (g/kg) 0.248 0.069 0.175 0.145 0.131 0.081 0.084 0.154 0.087 0.054 
TP (g/kg) 1.498"  1.390° 1.278° 0.9997 1.233  1.180° 1.042? 1.269 0.883*  0.872° 
SD (g/kg) 0.205 0.458 0.149 0.079 0.182 0.164 0.053 0.278 0.070 0.067 
TK (g/kg) 10.534? 9.961 8.231 7.243* 9.488? 8.036? 7.547} 8.480° 6.677"  6.201° 
SD (g/kg) 0.555 1.025 0.846 0.969 0.566 1.232 1.108 1.104 1.191 1.194 
C:N 16.793% 12.135" 17.577% 14.499% 21.895? 22.578" 15.880% 17.197 16.940 22.539° 
SD 1.012 1.467 1.307 2.574 1.858 2.775 1.992 2.033 2.089 0.175 
N:P 1.489" 0.590? 1.097% 0.853% 1.026% 0.888 0.867 0.854  0.972® 1.255% 
SD 0.254 0.083 0.195 0.135 0.083 0.212 0.130 0.167 0.223 0.240 
C:P 22.652" 14.248? 11.743? 25.019? 22.597? 15.380" 14.777 13.666° 13.573? 14.466° 
SD 4.322 2.647 2.007 5.502 2.817 2.208 1.991 2.163 2.105 2.549 


Note: SD, standard error. Different lowercase letters within a row indicate significant differences among different recovery time after 
fire and fire severities at P<0.05 level. LS, low-severity; MS, moderate-severity; HS, high-severity. 


The litter mass was significantly higher in the control (5.861 (+0.732) t/hm°) than those of the 


HS plots in 1974, 1983 and 1995 with the following order of control>LS>MS>HS (Fig. 1). The 
average litter masses for the LS, MS and HS plots of larch forest were 3.763, 3.303 and 2.512 t/hm2, 
respectively. For the plots in 1974, 1983 and 1995, the average litter mass decreased by 3.636, 
3.020 and 2.920 t/hm?, respectively (Fig. 1). 
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Fig. 1 Litter mass among the plots under different fire severities and recovery time after fire. Different 
lowercase letters mean significant differences among different recovery time after fire and fire severities at 
P<0.05 level. Error bars indicate standard errors. LS, low-severity; MS, moderate-severity; HS, high-severity. 


The SOC concentrations were positively related to the TBA (R?7=0.303; P<0.01). Significant 
positive relationships were observed between SOC and TBA (R?=0.270; P<0.01), litter mass and 
TBA (R7=0.270; P<0.01) (Fig. 2). 
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Fig. 2. Relationships between litter mass and total basal area (TBA; a), soil organic carbon (SOC) concentration 
and TBA (b), and SOC concentration and litter mass (c) 


3.2 Concentrations and contents of SOC, TN, TP and TK 


The concentrations of SOC, TN, TP and TK in the three soil layers (0—10, 10-25 and 25-70 cm) 
were highest in the control (Fig. 3). There were significant differences in SOC, TN, TP and TK 
concentrations between the HS plots and the control among the three soil layers. SOC, TN and TP 
concentrations in the MS plots were also significantly different from the control, except for the 
MS plot in 1983 in the 25-70 cm soil layer. However, the concentrations of SOC had significant 
differences in 1974 and 1983 under different fire severities only in the 0-10 cm soil layers (Fig. 
3a). TN concentrations in the 0-10 cm layer in 1974, the 10-25 cm layer in 1974 and 1983, and 
the 25-70 cm layer in 1974 under different fire severities showed similar trends (Figs. 3d—f). TP 
concentrations also showed similar trends, except for the plots in the 25-70 cm layer in 1995. The 
order of the concentrations of SOC, TN, TP and TK under different fire severities was as follows: 
control>LS>MS>HS. The order of concentrations of SOC, TN, TP and TK under different 
recovery time after fire was as follows: 1917>1974>1983>1995. 

In this study, SOC, TN and TP contents were higher in the control than in the other treatments, 
except for TK content (Fig. 4). The order of the SOC, TN, TP and TK contents under different fire 
severities was as follows: control>LS>MS>HS, except for the TK contents in the 25-70 cm layer. 
The order of SOC and TN contents under different recovery time after fire was as follows: 
1917>1974>1983>1995. TP and TK contents showed no significant differences under different 
recovery time after fire and fire severities except for the 0-10 cm soil layer in 1993 and the 10-25 
cm layer in 1974. 


3.3 SEM analysis 


The SEM fitting was used to analyse the effect of recovery time after fire and fire severity on 
SOC and soil nutrients. Average concentrations of SOC, TN, TP and TK were used for analysis, 
and the control data were excluded (P=0.121; Fig. 5). Results suggested that recovery time after 
fire and fire severity have direct impacts on SOC and soil nutrients. Increasing fire severity 
appeared to decrease TN (—0.34), SOC (—0.24), TK (—0.23) and TP (—0.21) concentrations. 
Positive effects of recovery time after fire on TN, SOC, TK and TP concentrations were observed, 
and the standardized path coefficients were 0.49, 0.29, 0.25 and 0.18, respectively. 


4 Discussion 


Recovery time after fire and fire severity exerted significant effects on stand structure of larch 
forest and soil physical and chemical properties in the KNNR, Northwest China. The control 
generally had higher TBA and litter mass of larch forest and had higher concentrations in SOC, 
TN, TP and TK, compared with burned stands. Among stands that had been disturbed within the 
last 43 a, increasing recovery time after fire and decreasing fire severity generally led to increased 
TBA, litter mass, SOC, TN, TP and TK concentrations; however, not all the differences were 
generally significant. 
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Fig. 3 Soil organic carbon (SOC), total nitrogen (TN), total phosphorus (TN) and total potassium (TK) 
concentrations in different soil layers under different fire severities and recovery time after fire. LS, low-severity; 
MS, moderate-severity; HS, high-severity. Different lowercase letters show significant differences among 
different recovery time after fire and fire severities at P<0.05 level. Bars indicate standard error. 


4.1 Soil property and stand structure 


The significantly different C:N, N:P and C:P ratios revealed variations of SOC, TN and TP 
between post-fire and pre-fire conditions. Roscoe et al. (2000) and Francos et al. (2018) also 
identified significant differences in the soil C:N ratio several decades after wildfire. Fire severity 
affects stand structure of forest and soil recovery, and its effects may last many years after a fire 
(Dzwonko et al., 2015; Miquelajauregui et al., 2016). In this study, litter mass increased with 
increasing recovery time after fire, and the linear regressions between litter mass and SOC, TBA 
and SOC, and litter mass and SOC were positive, which may indicate that the input of organic 
matter is less than litter production and decomposition (Hu et al., 2013). These findings suggested 
that fire severity had long-term impacts on litter mass and TBA and TBA impacted SOC through 
litter mass losses. The litter mass of the LS plots in 1974 was 4.710 t/hm?, which was less than 
that of the control (5.861 t/hm?). Fire eliminated most of the litter, and the losses caused by HS 
fire were greater than those of LS and MS ones, and the recovery rate of litter in the LS plot was 
the highest. There was less litter mass in the plots with HS fire than in the plots with MS and LS 
fires, which is consistent with the findings of Pausas et al. (2002). It took more than 43 a for 
forest litter to recover to the pre-fire level, depending on forest structure and fire severity. The 
relationship between litter mass and TBA was positive, indicating that leaf litter production was 
lower in the youngest plots (Quintero-Gradilla et al., 2015). TBA decreased, but stand density 
increased with increasing fire severity, supporting the findings that fire strongly alter the forest 
structure (Kasischke and Johnstone, 2005; Laughlin and Fulé, 2008; Martin et al., 2012). 
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Fig. 4 Soil organic carbon (SOC), total nitrogen (TN), total phosphorus (TP) and total potassium (TK) contents 
in different soil layers under different fire severities and recovery time after fire. LS, low-severity; MS, 


moderate-severity; HS, high-severity. Different lowercase letters show significant differences among different 
recovery time after fire and fire severities at P<0.05 level. Bars indicate standard errors. 


0.09 


Fig. 5 Result of structural equation model (SEM). Time, recovery time after fire; Severity, fire severity; SOC, 
soil organic carbon; TN, total nitrogen; TP, total phosphorus; TK, total potassium. 


4.2 Concentrations and contents of SOC and soil nutrients 


SOC and TN concentrations were found to increase with increasing recovery time after fire, 
which is consistent with previous studies, indicating that soils need time to recover after a fire 
(Simansky, 2015; Francos et al., 2018). It was reported that soil C and N are reduced due to 
volatilization, charring and oxidation under HS fire (Giovannini et al., 1988; Rovira et al., 2012). 
We found that SOC, TN, TP and TK concentrations in the 0-10 and 10-25 cm soil layers need 
more than 43 a to recover to the pre-fire levels. The variations in TP and TK concentrations in 
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different soil layers under different recovery time after fire were not consistent, likely because soil 
properties can be influenced by recovery time after fire and soil depth (Verma and Jayakumar, 
2018). In the 0-10 cm layer, the SOC and TN concentrations in the plots with LS fire were higher 
than those in the plots with MS and HS fires, and the corresponding concentrations in the control 
were higher compared with those under shorter recovery time after fire. After a fire, SOC 
concentration significantly decreased (Bell and Binkley, 1989; Murphy et al., 2006; Badia et al., 
2012; Alcaniz et al., 2016), and TN concentration could also decrease. A wildfire at Donnelly 
Flats, Alaska, USA resulted in the losses of up to 50% of soil C and 33% of soil N stocks (Neff et 
al., 2005). The results of the present study are consistent with those of Debano (2000), who found 
that the strongest responses in SOC and TN concentrations existed in the 0-10 cm soil layer. The 
effects of fire on SOC, TN, TP and TK concentrations decreased gradually with increasing soil 
depth, demonstrating that the impacts of fire severity were weaker in deeper soil than in surface 
soil (Fernandez et al., 1997; Certini, 2005; Wang et al., 2013). 

We observed that SOC contents were significantly different in the 0—10 cm soil layer for the 
same year of fire occurrence, suggesting that fire severity greatly impacted SOC contents. The 
SOC contents in the plots with longer recovery time after fire were higher than those of with 
shorter recovery time after fire. Fernandez et al. (1997) also found that more than 50% of the 
organic matter in the surface soil humus was lost after a wildfire. In the moderate- and long-terms, 
fire-affected areas have significantly lower total carbon content than the control plot (Francos et 
al., 2018). Johnson et al. (2005) and Kaye et al. (2010) observed that total carbon content failed to 
recover to the pre-fire levels for 10-20 a after a wildfire. 

The results of the present study suggested significant differences in TN contents among the 
plots with different fire severities and the same recovery time after fire in the 0-10 and 10-25 cm 
soil layers. This indicated that HS fire had more significant impacts on TN contents. The TN 
contents in the control plot were higher than those of the plots with different recovery time after 
fire, suggesting that it took a long time to recover to the pre-fire TN levels (Moghaddas and 
Stephens, 2007). The decrease in soil TN contents with increasing burning temperature may be 
explained by the reduction in organic N contained in organic matter (Huang and Boerner, 2007) 
and direct volatilization (Burgoyne and Deluca, 2009). The significant differences among TN 
contents in three soil layers with HS fire and the control indicated that TN content was sensitive 
to fire severity. 

The TP contents in the 0-10 and 10-25 cm soil layers decreased with increasing fire severity, 
which is in line with the results of Certini (2005). Duguy et al. (2007) observed a decrease in total 
and available P contents in soil after a fire in eastern Spain. Moreover, direct volatilization and 
oxidation (Raison, 1985) and nutrient loss as particulates in smoke may cause losses of TP 
content (Certini, 2005). We did not observe any significant variation in TK contents among the 
burned plots, which is different from the Mehdi et al. (2012). The trend of TK content might be 
attributed to the similar value of bulk density among the burned plots. 


4.3 Relationships of fire disturbance factors with concentrations of SOC and soil nutrients 


The SEM was applied for analyzing the relationships of fire disturbance factors and 
concentrations of SOC and soil nutrients. The results suggested that SOC, TN, TP and TK 
concentrations increased with decreasing fire severity and increasing recovery time after fire. 
Forest floors are a major reservoir of soil N that is substantially decreased under a fire by the 
burning of litter, duff and soil (Driscoll et al., 1999; Simard et al., 2001). Giovannini et al. (1988) 
observed that fire caused total destruction to partial scorching of organic matter, including 
volatilization, charring and oxidation. Verma and Jayakumar (2018) found that losses of 
phosphorus by volatilization or leaching were small. The retrogradation of potassium was 
attributable to precipitation and adsorption after a fire and then increased by the colonization of 
the burned areas by vegetation during recovery (Alcaniz et al., 2016). The above reports agreed 
with our results that the sensitivity order of SOC and soil nutrients to fire severity and recovery 
time after fire was TN, SOC, TK and TP, indicating that fire severity and recovery time after fire 
play important roles in soil recovery in burned areas. 


5 Conclusions 


Fire severity and recovery time after fire affected the stand structure, litter mass of larch forest 
and soil properties in the KNNR, Northwest China. Stand structure, litter mass, SOC and soil 
nutrients in the 0-10 cm and 10-25 cm soil layers need more than 43 a to recover to the pre-fire 
levels, and the recovery rate varies with fire severity and recovery time after fire, with HS fire 
needing the longest time to recover. The slow recovery rate may result from the low metabolic 
rate of the boreal forest ecosystem. Fire severity plays a crucial role in the carbon and nutrition 
cycles of the boreal forest ecosystem, and it should be noted that fire severity will likely increase 
in the future due to climate change. Therefore, future fire in this area should be controlled within 
the range of low-severity. It is recommended that field surveys should be conducted immediately 
at least one to three years after fire to analyze the short-term effects of fire in the KNNR. 
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